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ABSTRACT 

We present a compilation of high time-resolution photometric observations of quiescent 
soft X-ray transients obtained with the acquisition camera of Gemini-South. A0620- 
00 was observed with a short cycle time and high precision. Superimposed on the 
ellipsoidal modulation we find several prominent flares together with weaker continual 
variability. The flares seen sample shorter timescale than those reported in previous 
observations, with rise times as low as 30 s or less; most flares show unresolved peaks. 
The power density spectrum (PDS) of A0620-00 appears to exhibit band-limited noise 
closely resembling the X-ray PDS of black hole candidates in their low states, but 
with the low-frequency break at a lower frequency. X-ray Nova Muscae 1991 shows 
much larger amplitude flares than A0620-00 and if a break is present it is at a lower 
frequency. XN Vel 1993 shows very little flaring and is, like A0620-00, dominated 
by the ellipsoidal modulation. We discuss the possible origins for the flares. They 
are clearly associated with the accretion flow rather than an active companion, but 
whether they originate in the outer disc, or are driven by events in the inner region is 
not yet resolved. The similarities of the PDS to those of low/hard state sources would 
support the latter interpretation, and the low break frequency is as would be expected 
if this frequency approximately scales with the size of an inner evaporated region. We 
also report the discovery of a new variable star only Idarcsec from XN Mus 1991. 
This appears to be a W UMa star, with an orbital period of about 6hrs. 

Key words: accretion, accretion discs - binaries: close - stars: individual: V616 Mon, 
GU Mus, MM Vel 


1 INTRODUCTION 

Soft X-ray transients (SXTs), also referred to as X-ray novae 
and black hole X-ray transients, are low-mass X-ray binaries 
in which long periods of quiescence, typically decades, are 
punctuated by very dramatic X-ray and optical outbursts, 
often accompanied by radio activity (Tanaka & Shibazaki 
1996; Cherepashchuk 2000). In outburst a number of X-ray 
spectral states are seen, most commonly the high/soft state 
and the low/hard state. An intermediate and very high state 
have also been identified. In the high/soft state. X-ray emis¬ 
sion is dominated by thermal emission from an accretion 
disc extending to close to the last stable orbit around a 
black hole. In the low/hard state, the inner disc is believed 
to be truncated and emission appears to originate from an 
extended corona. Direct support for this picture is provided 
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by the low/hard state source, XTE J1118-I-480 (Hynes et al. 
2000), in which an inner disc radius of at least SORsch, and 
probably ^ 350Rsch, is required (McClintock et al. 2001; 
Chaty et al. 2002). A similar scenario is advanced by vari¬ 
ous advective models for the quiescent state (see Narayan, 
Garcia & McClintock 2001 and references therein), but with 
the disc truncated at larger radii, 10^-10® Rsch- Attempts 
have been made to unify these spectral states within the ad¬ 
vective picture (Esin et al. 1997), and Esin et al. (2001) did 
achieve some success in fitting the broad band spectrum of 
the low/hard state source XTE J1118-1-480 with an advec¬ 
tive model. 

The states of SXTs are classified by their X-ray timing 
properties as well as by their spectra. With the exception 
of the quiescent state, these have been well studied (e.g. 
van der Klis 1995; Wijnands & van der Klis 1999). The 
high/soft state shows a low level of red noise, with no de¬ 
tected low frequency break. The low/hard state and very 
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high state exhibit a higher level of band-limited noise, with 
a low-frequency break at ~0.02-30 Hz, and sometimes su¬ 
perposed QPOs. Band-limited noise also appears to be seen 
in the Seyfert 1 galaxy NGC 3516 (Edelson & Nandra 1999), 
but with a much lower cutoff, 4 x 10“^ Hz. This suggests an 
approximate scaling with black hole mass, and presumably 
with the scale of the accretion region. The AGN data do 
not extend to low enough frequencies to be confident that 
the break is analogous to the low frequency break in the 
low/hard state, however; it could actually correspond to a 
higher frequency turnover (Uttley 2002, priv. comm.). 

Similar properties might be expected for quiescent 
SXTs, as the structure of the flow is believed to be sim¬ 
ilar to that in the low/hard state. Observations of quies¬ 
cent state variability are much more difficult, however. Sub¬ 
orbital variability is known to be present at X-ray energies 
in the brightest source, V404 Cyg (Wagner et al. 1994; Kong 
et al. 2002), but even this is faint; only ~ 0.15photonss“^ 
are detected with Chandra. Variability can be more effec¬ 
tively studied in the optical where reasonable count rates 
are possible, and several photometric (Haswell et al. 1992; 
Pavlenko et al. 1996; Zurita, Gasares & Shahbaz 2002a) and 
spectrophotometric (Hynes et al. 2002) studies have been 
performed. The origin of the variability remains uncertain, 
however, with plausible possibilities including direct optical 
emission from an advective region, reprocessed X-ray vari¬ 
ability, magnetic reconnection events in the disc, and flick¬ 
ering from the accretion stream impact point. It may be 
that a combination of these factors are important on differ¬ 
ent timescales, and that not all objects are dominated by 
the same source of variability. It is therefore important to 
perform a comparative study of the class as a whole to de¬ 
termine if there is just one type of variability or two or more 
with distinct characteristics. By doing this we can hope to 
isolate the contribution, if any, from variations in the inner 
flow and hence probe its nature. 

Even in the optical, count rates from existing data are 
quite low, limiting the time-resolution achieved, and severely 
compromising data quality for fainter objects. As we are 
studying aperiodic variability, the only good solution to this 
problem is to increase the count rate with a larger aperture 
telescope. Consequently we have embarked upon a survey 
of fast variability with the Acquisition Camera on Gemini- 
South. Our main goals are to explore the bright targets on 
faster timescales than previously possible, and to study vari¬ 
ability effectively even in faint objects. We present here high 
time resolution data on the bright prototypical SXT, A0620- 
00 (=V616 Mon), together with lower resolution lightcurves 
of the fainter objects X-ray Nova Muscae 1991 (GU Mus), 
and X-ray Nova Velorum 1993 (MM Vel). 


2 OBSERVATIONS 

Photometric observations of the three quiescent SXTs were 
obtained in service mode with the Acquisition Camera (Acq- 
Cam) on Gemini-South on 2001 December 15 and 2002 Jan¬ 
uary 11 and 15. A V filter was used with exposures between 6 
and 60 s. The minimum dead-time was primarily dictated by 
the CCD readout time, 1.7 s, but a variable additional delay 
was involved in data transfer. The typical actual dead time 
between exposures was about 2.2 s. Full details are given in 


Table Routine bias and dark subtraction and flat-fielding 
corrections were applied to all the data before distribution 
and appear satisfactory. 

A0620-00 and XN Mus 1991 were relatively bright com¬ 
pared to the background and are not crowded (no significant 
contaminating stars closer than 5arcsec.) These were there¬ 
fore straightforward to analyse. Lightcurves of each object 
and several comparison stars were extracted using aperture 
photometry. For A0620-00, a weighted mean of the magni¬ 
tudes of two nearby brighter comparisons was used for dif¬ 
ferential photometry; for XN Mus 1991 four brighter com¬ 
parisons were used. In each case, two fainter nearby non¬ 
variable comparisons of the same brightness as the target 
were used to check the photometric accuracy. The photo¬ 
metric apertures (0.7 arcsec for A0620-00, and 0.6 arcsec for 
XN Mus 1991) were chosen to minimise the variance in the 
lightcurve of one of the faint comparison stars. From these 
comparisons we estimate la accuracies of 0.8percent per 
exposure for the flux lightcurve of A0620-00 and 1.7 percent 
for XN Mus 1991. 

XN Vel 1993 was more problematic as it lies close 
to a brighter star. Fortunately, the seeing was good; 0.5- 
0.7 arcsec FWHM for most of the run with a few images at 
the beginning with FWHM up to 0.9 arcsec. An additional 
complication was that the immediate field of XN Vel 1993 
was located close to the seam between the two halves of the 
detector. To minimise the effect of this on sky estimates, we 
interpolated over the gap. XN Vel 1993 itself was not affected 
and we did not use any stars that were on the seam either as 
PSF or comparison stars. Photometry of XN Vel 1993 and 
a number of comparison stars was obtained by PSF fitting 
using the iraf|^ implementation of daophot. PSF fitting 
with AcqGam data is difficult, as the camera optics yield a 
position dependant PSF. This can largely be dealt with by 
choosing many PSF stars surrounding the objects of interest 
and modelling position dependence. We restricted the anal¬ 
ysis to the immediate region of XN Vel 1993, with 7 suitable 
PSF stars, and allowed only linear variations (with respect 
to X and y) in the PSF; there were too few PSF stars to use 
a higher order model. This process did leave small residuals 
in the core of the PSF of the brighter stars, but the wings 
are well subtracted. Fortunately, in most of our images, XN 
Vel 1993 is sufficiently resolved from the brighter star that 
the latter should not be a significant problem. For the im¬ 
ages at the beginning of the run, however, the poorer seeing 
caused problems and these images were excluded from the 
extracted lightcurve. Excluding these points, the extracted 
lightcurves were insensitive to whether a variable PSF was 
used or not, and to the adopted fitting radius. Consequently 
we believe that the lightcurve of XN Vel 1993 is not signifi¬ 
cantly contaminated by difficulties in the fitting. Differential 
photometry was performed with respect to six nearby com¬ 
parison stars, all brighter than XN Vel 1993 and relatively 
isolated. Another comparison star of similar brightness to 
XN Vel 1993 was extracted. Both the formal errors and the 
scatter in the comparison lightcurve give ala uncertainty 
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Table 1. Log of Gemini-South AcqCam observations. 


Object 

Date 

UT range 

Exposures 

Airmass 

Seeing 

(V) 

A0620-00 

2001 Dec 15 

02:26-06:30 

1705 X 6s 

1.2-1.7 

0.6-1.0 

18.1 

XN Mus 1991 

2002 Jan 11 

04:32-08:35 

219 X 60s 

1.3-1.6 

0.5-1.1 

20.3 

XN Vel 1993 

2002 Jan 15 

04:23-08:15 

425 X 30 s 

1.0-1.3 

0.5-0.7 

21.7 



Orbital phase 

Figure 1. Lightcurve of A0620—00. The abscissa is orbital phase 
with respect to the ephemeris of Leibowitz et al. (1998). A small, 
~ 0.05, phase offset is present. The comparison star has the same 
average brightness as A0620—00, but has been offset downwards 
by 0.15 units for clarity. The fits to the ellipsoidal modulation are 
based on a double sinusoid (fundamental on the orbital period 
plus first harmonic). In the equal maxima case, the phases of 
the sinusoids have been fixed to produce a lightcurve like a pure 
ellipsoidal case. In the unequal case, the phasing is allowed to 
float to better fit the data, as appropriate if the orbital lightcurve 
is distorted by other effects. Letters indicate the regions expanded 
in Fig. |. 


of 0.9 percent for this star. Since XN Vel 1993 is blended, 
it is subject to larger errors; the formal error estimate is 

1.1 percent. 

For all three objects, an approximate photometric cal¬ 
ibration was applied with respect to several standard fields 
observed at low airmass on the same night. Colour terms 
were neglected as the targets were observed in a single band 
only. 

3 LIGHTCURVES 

3.1 A0620-00 

The lightcurve of A0620-00 is shown in Fig. It is clearly 
dominated by the ellipsoidal modulation due to the distor¬ 
tion of the companion star (c.f. McClintock & Remillard 
1986). Since this is not associated with accretion variability 
we should remove this to isolate the flares. We do this in 
a comparable way to Zurita et al. (2002a), by approximat¬ 
ing it with two sine waves at the orbital frequency and its 


hrst harmonic. Most of the modulation is at the hrst har¬ 
monic but the fundamental provides for a variation in the 
minima and/or maxima. For an ideal pure ellipsoidal modu¬ 
lation the maxima are equal, but the minima differ. In prac¬ 
tise, lightcurves of quiescent SXTs in general, and A0620- 
00 in particular, often exhibit unequal maxima as well (e.g. 
Haswell 1996; Leibowitz, Hemar & Orio 1998). This may 
be due to distortion of the orbital lightcurve by light from 
the stream impact point, starspots, and/or persistent super¬ 
humps; the latter definitely appear to be seen in one source, 
XTE J1118-I-480, albeit in the last stages of outburst de¬ 
cline rather than true quiescence (Zurita et al. 2002b). We 
consider both cases; for the pure ellipsoidal lightcurve we 
fix the relative phases of the sine waves to produce equal 
maxima. We also allow the phases to vary independently. 
We ht using an iterative rejection scheme to approximately 
ht the lower envelope of the lightcurve, by rejecting points 
more than 2 a above the fit, then refitting. This is repeated 
until no new points are rejected. We show the results of the 
fitting in Fig. Allowing the maxima to differ does improve 
the fit somewhat, but the difference is not dramatic. 

Superposed on the ellipsoidal modulation are many 
rapid flares; the strongest are shown in Fig. Similar flares 
in A0620-00 have been reported by Haswell (1992) and Zu¬ 
rita et al. (2002a). We sample events of shorter duration, 
however; both of the previous studies had a time-resolution 
of 30-40 s, which would barely have resolved our shortest 
events. For the most prominent and distinct flares shown 
in Fig. ^ we have estimated some characteristics of each 
flare; the peak amplitude, equivalent duration (c.f. Zurita 
et al. 2002a), and rise and decay e-folding timescales, based 
on an exponential ht. These are not a representative sam¬ 
ple, as they are selected to be the most extreme, best dehned 
events. It can be seen that these flares lie at the extreme-low 
end of the distribution of equivalent durations presented by 
Zurita et al. (2002a). E-folding timescales are typically 30- 
80 s, although one event (hare 5) rises much more rapidly 
than this. There is no consistent asymmetry to the hares, 
although individual events may be asymmetric; again the 
most extreme behaviour is shown by hare 5. We have also 
characterised the activity level of the lightcurve using similar 
nomenclature to Zurita et al. (2002a) in Table p. 


3.2 X-ray Nova Mus 1991 

The lightcurve of XN Mus 1991 is shown in Fig. Large am¬ 
plitude aperiodic variability is dominant, to the extent that 
any contribution from an ellipsoidal variation is not obvi¬ 
ous in these data, and cannot be htted. Ellipsoidal modula¬ 
tions have previously been reported with a full amplitude in 
B V of ~0.2-0.35 mag (Remillard, McClintock & Bailyn 
1992; Orosz et al. 1996). Their apparent absence is probably 
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Table 3. Properties of Gemini lightcurves. the spectroscopic veiling, is the contribution due to the non-variable disc light, "zf 

is the mean flare flux and cr^ its standard deviation, both expressed as a fraction of the mean flux, a* = azjv'^ and rj is the fraction of 
the average veiling due to the flares. 


Object 

'^obs 



O'z 


7 

A0620-00 

6±3 (Ha), 17 ±3 (H/3) 

12 ±3 

0.009 

0.015 

0.12 

0.07 

XN Mus 1991 

~ 54 (5000 A), 15 ± 6 (6400 A) 

26 ±4 

0.23 

0.14 

0.56 

0.47 

XN Vel 1993 

65 ± 5 (6300 A) 

64 ±5 

0.024 

0.023 

0.04 

0.04 


Table 2. Properties of selected individual flares from A0620—00. 
The peak is defined as a fraction of the mean level. The equiva¬ 
lent duration is the total counts in the flare divided by the mean 
counts. The rise and decay times are e-folding times from expo¬ 
nential fits to the rising and decaying segments of the flare profile. 


Flare 

Peak 

Equivalent 

Rise 

Decay 

No. 


duration (s) 

time (s) 

time (s) 

1 

0.05 

2.8 

30 

76 

2 

0.10 

9.5 

29 

69 

3 

0.05 

4.0 

80 

39 

4 

0.06 

6.2 

7 

81 

5 

0.12 

7.7 

34 

32 


due to a combination of large flare amplitude, small ellip¬ 
soidal amplitude and relatively long orbital period. We show 
in Fig. ^ the expected modulation, assuming the ephemeris 
of Shahbaz et al. (1997) and a full amplitude of 0.27mag. 
It clearly is not consistent with the data. However, there 
may be an error in extrapolating the ephemerides over sev¬ 
eral years, as appeared to be the case for A0620-00; using 
the formal error estimate of Shahbaz et al. (1997) this cor¬ 
responds to an uncertainty of 0.06 in phase at our epoch. 
Fig. I also shows the effect of including a phase offset of 0.12 
(i.e. 2 a), which is not unreasonable; the agreement with the 
lower envelope of the data is now acceptable, and with this 
allowed for our data are consistent with an ellipsoidal am¬ 
plitude comparable to that previously reported. There are 
too few points defining the lower envelope to reliably fit a 
model with both the phase and amplitude variable, however, 
so we neglect the relatively small ellipsoidal contribution in 
subsequent sections unless explicitly noted. Since the flaring 
in XN Mus 1991 is of such large amplitude, relative to the 
ellipsoidal contribution, this should not introduce a large 
error. 

Since the time-resolution was much less than for A0620- 
00, we cannot study rise and decay timescales of flares 
so readily although relatively fast timescales are clearly 
present. For example the large and extended flare at the 
end of the observation represents a net increase of about 
40 percent, with an e-folding time of ^ 50 s. As for A0620- 
00, the activity level of the lightcurve is characterised in 
Table ^ to allow comparison with the results of Zurita et al. 
(2002a). 



0 500 1000 0 500 1000 

Elapsed time (s) 


Figure 2. Close-up view of a selection of flares in A0620—00. 
The ellipsoidal modulation has been subtracted and replaced by 
its average value. Time and flux scales are the same for easy 
comparison of flare amplitudes and durations. The zero points of 
elapsed time are arbitrary. Numbered flares have their properties 
summarised in Table H. 


3.3 X-ray Nova Vel 1993 

The lightcurve of XN Vel 1993 is shown in Fig. ^ phased on 
the new ephemeris of Gelino (priv. comm.). This is domi¬ 
nated by an apparent ellipsoidal modulation (c.f. Shahbaz et 
al. 1996). Superposed on the ellipsoidal modulation do ap¬ 
pear to be some flares. The flare amplitude is much less than 
in XN Mus 1991, but is comparable to that seen in A0620- 
00. Parameters of the detrended variability are summarised 
in Table ^ 

4 POWER DENSITY SPECTRA 
4.1 A0620-00 

To quantify the range of timescales present, we calculate 
a power density spectrum (PDS), after removing the fit¬ 
ted ellipsoidal modulation from the lightcurve. Since the 
sampling was not perfectly uniform we calculate a Lomb- 
Scargle periodogram (see Press et al. 1992 and references 
therein) and normalise it in the same way as is common 
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Orbital phase 


Figure 3. Lightcurve of XN Mus 1991. The abscissa is orbital 
phase with respect to the ephemeris given in Shahbaz, Naylor & 
Charles (1997). The comparison star has the same average bright¬ 
ness as XN Mus 1991, but has been offset downwards by 0.3 units 
for clarity. The dashed line indicates the expected ellipsoidal mod¬ 
ulation. This is not consistent with the data, but allowing an 0.12 
phase offset (solid line) improves the situation. 



Frequency (Hz) 


Figure 5. PDS of A0620-00 and XN Mus 1991. The white noise 
component, fitted at the highest frequencies, has been subtracted. 
The fitted red noise slopes are —1.52 for A0620—00 and —1.56 
for XN Mus 1991. The break frequency in A0620-00 is at 9.5 X 
10““* Hz, corresponding to a timescale of ~ 20 min. For XN Mus 
1991, filled squares indicate the PDS of the raw lightcurve, the 
open squares indicate the effect of removed the representative 
ellipsoidal modulation shown in Fig. M 



Orbital phase 


Figure 4. Lightcurve of XN Vel 1993. The abscissa is orbital 
phase with respect to the updated ephemeris derived by Gelino 
(priv. comm.). The comparison star has almost the same average 
brightness as XN Vel 1993, but has been offset downwards by 0.3 
units for clarity. The fits to the ellipsoidal modulation are based 
on a double sinusoid model exactly as described for A0620-00. 


for Fourier transform PDS. Following the suggestion of Pa- 
padakis & Lawrence (1993), we bin and fit the PDS in loga¬ 
rithmic space, i.e. each bin is evaluated as (logp) rather than 
log (p) and we then fit to the values of (logp). Errors on the 
binned logarithmic power are estimated from the standard 
deviation of points within the bin. White noise has been 
estimated and subtracted by fitting the highest frequencies 


with a white noise plus red noise model, but the white noise 
was not large, since the photometric precision was good. 

The derived PDS is shown in Fig. ^ At higher fre¬ 
quencies we see a well defined power-law. Below 1 mHz 
(i.e. (20min)“^), this appears to flatten. The overall bro¬ 
ken power-law form is strikingly similar to the X-ray PDS 
of SXTs in the low/hard state, as will be discussed in Sec¬ 
tion 3.4. To characterise it numerically, we fitted it with a 
model which has a red noise power-law above a break fre¬ 
quency and is flat below that. We derive a power-law slope 
of —1.52, similar to low/hard state SXTs and a break fre¬ 
quency at 9.5 X 10“^^ Hz. 

There are a number of potential pitfalls with this analy¬ 
sis. In detrending the ellipsoidal modulation we may have re¬ 
moved low frequency accretion variability as well and hence 
flattened the PDS artificially. Also, with only 4hrs of data 
the shape of the low-frequency PDS will be subject to fluc¬ 
tuations due to the individual realisations of the spectrum. 
Finally there may be some aliasing problems at high frequen¬ 
cies due to the deadtime between exposures and the non¬ 
uniformity of the sampling. To attempt to quantify these 
uncertainties we create simulated lightcurves with exactly 
the same sampling and integration times as the real data. 
We begin with a modelled ellipsoidal lightcurve, to ensure 
that the double-sine approximation used does not provide 
an exact match and that any systematic error introduced by 
the inadequacy of this model is reproduced by the simula¬ 
tion. This was calculated assuming Mx = 10 Mq, q = 0.067, 
i = 40°, iodise = 0.45 Rli. These are arbitrary choices and 
are only intended to produce a representative simulated 
lightcurve; this model was not used to detrend the data. 
With these parameters a reasonable agreement was obtained 
with the observations, however. To this was added a model 
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noise lightcurve, calculated using the method of Timmer & 
Konig (1995). Several models for the PDS were tested, bro¬ 
ken power-laws with a break close to that observed, 10~® Hz, 
at the edge of the useful coverage, 10“^^ Hz and well out¬ 
side the range sampled, 10“® Hz. In each case the power-law 
part of the spectrum was matched to the observed slope 
and normalisation. Lightcurves were created with 1 s time 
resolution and binned up so that high frequency variabil¬ 
ity would be aliased correctly. We included lower frequency 
variations than are well sampled by our observations so that 
red noise leaks would have an effect. Finally Gaussian white 
noise was added using the errors derived from photometry; 
for the comparison star these do represent the scatter in the 
data accurately so the noise level should be comparable to 
the real data. 

For each model of the PDS we calculate 1000 simulated 
lightcurves and analyse them in exactly the same way as 
we did the real data. We create individual PDS with the 
same logarithmic frequency binning used for the data, and 
then compare the range of PDS obtained with the obser¬ 
vations. We find that if a broken power-law model fitting 
the observations is input, then the simulated PDS do match 
the observations (Fig. ^), i.e. we would not expect such a 
PDS to be substantially distorted from the input form. If 
the break is dropped to 10“^ Hz then the observed range of 
the PDS corresponds to an unbroken power-law. The out¬ 
put of the simulations preserves this form fairly well, al¬ 
though there is some flattening at low frequencies due to 
removal of some variability at the orbital frequency and its 
first harmonic. This flattening is at lower frequencies than 
the observed break in the PDS of A0620-00, and this model 
does not agree with the data (Fig. ^). Only 1 simulation 
out of 1000 produced sufficient flattening that all bins below 
10~® Hz have a power density below 0.1 Hz“^ as observed. 
We also ran simulations with even lower frequency breaks, at 
10“® Hz, to allow for the possibility of distortion by strong 
red noise leaks. These also provided a very poor agreement 
with the observations. We can therefore reject a model in 
which the intrinsic PDS breaks at or below the edge of the 
observed range, and conclude that the break in the PDS is 
probably real. 

4.2 X-ray Nova Mus 1991 

We calculated a PDS for XN Mus 1991 in exactly the same 
way as for A0620-00. It is also shown on Fig. The PDS of 
the raw lightcurve shows no sign of the break seen in A0620- 
00, instead possessing an unbroken power-law PDS. The 
slope is very similar to that of A0620-00, however, —1.56. 

Of course, just as subtracting a fitted ellipsoidal mod¬ 
ulation from A0620-00 may distort the PDS if some aperi¬ 
odic power is removed as well, so also an uncorrected ellip¬ 
soidal modulation could distort the PDS of XN Mus 1991, 
although in the opposite sense. To test for this, we also tried 
subtracting the representative (not fitted) ellipsoidal modu¬ 
lation shown in Fig. ^from the lightcurve before calculating 
the PDS. Since we cannot be confident that this is an ac¬ 
curate model, the resulting PDS may be no more correct 
that without subtraction, but the differences between them 
should indicate the sensitivity to the ellipsoidal contribution. 
We see that a break does emerge, but at a lower frequency 
than seen in A0620-00. Without being sure of the ellipsoidal 



10 “® 10 “^ 10 “' 
Frequency (Hz) 


Figure 6. Simulated PDS of A0620-00 compared with the data. 
The solid line in each panel indicates the average of 1000 simu¬ 
lations. The dashed lines indicate 1 a confidence regions for indi¬ 
vidual points, a) Our best fit broken power-law model, indicating 
that this is consistent with the data, and systematics due to sub¬ 
traction of the ellipsoidal modulation and red noise leaks do not 
distort it substantially from the input model, b) A case with sim¬ 
ilar slope, but the break moved to the edge of the observable 
range, 10“^ Hz. This is clearly a poor fit. 


contribution we cannot claim that this break is real, but we 
can say that any break in the PDS of XN Mus 1991 must 
be at < 3 X 10“"^ Hz. Removing the ellipsoidal contribution 
also flattens the PDS somewhat, as it otherwise introduces 
a red leak across a range of frequencies which steepens the 
PDS. 


4.3 X-ray Nova Vel 1993 

We attempted to calculate a PDS for XN Vel 1993, but the 
combination of low count rate and low variability amplitude 
meant that it was not well-defined over a useful frequency 
range. There is little prospect of substantially improving on 
this, unless XN Vel 1993 were caught in a state showing a 
higher level of flaring. 
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5 A SERENDIPITOUS VARIABLE STAR 
DISCOVERY 

One of the potential comparison stars we examined for XN 
Mus 1991 was revealed to be itself a variable. This is about 
Idarcsec SSW of XN Mus 1991, at RA ll‘'26™25!3 ± (f.2, 
Dec —68°40'44'.'5 ± iVO (J2000), based on interpolation be¬ 
tween 18 surrounding USNO A2.0 stars (Monet et al. 1998). 
From an approximate photometric calibration we estimate 
ly) = 20.5, where the average is obviously over the observed 
lightcurve, not a whole cycle. It displays a smooth modula¬ 
tion rather than erratic variability, reminiscent of W UMa 
contact binaries (Fig. From the maximum and minimum 
which are fully observed, a full amplitude of 0.37 mag is mea¬ 
sured. The later incomplete maximum, however, suggests an 
asymmetric higher peak, and a full amplitude of > 0.40 mag. 
Of course, if the other minimum is deeper then the full am¬ 
plitude will be even larger. The period of the modulation is 
either 2.87±0.02 hrs if single peaked, or 5.73±0.03 hrs if dou¬ 
ble peaked. The latter seems more likely, given the unequal 
maxima observed, and is a typical period for a W UMa star 
(Maceroni & van’t Veer 1996). The asymmetric maxima are 
quite common to this class of objects (e.g. Davidge & Milone 
1984 and references therein) and are often attributed to star 
spots on these active late type binaries. The difference we 
see, > 0.037 mag, is comparable to that seen in other sys¬ 
tems, and indeed the same difference as in the prototype, 
W UMa (e.g. Maceroni & van’t Veer 1996). In short, it is 
most likely that this variable is a W UMa binary, or a re¬ 
lated type. A more comprehensive study, based on several 
cycles of variability, would be needed to draw a more confi¬ 
dent conclusion. Such a study would likely be possible with 
existing images already obtained to study ellipsoidal mod¬ 
ulations in XN Mus 1991, as most such images should also 
include the new variable. 


6 DISCUSSION 

6.1 What determines the variability amplitude? 

Zurita et al. (2002a) compared observed levels of variabil¬ 
ity with various system parameters. The only correlation 
they found was with the binary inclination, but as they ex¬ 
plain, this is somewhat misleading. This is because it ac¬ 
tually arises from a more fundamental correlation with the 
veiling; the fractional variability correlates with the frac¬ 
tional contribution of disc light to the total. This is a very 
sensible result, but not trivial, since it indicates that the 
variability is associated with the disc, as expected, and not 
with activity on the companion star. 

We illustrate the correlation using both our data and 
that of Zurita et al. (2002a) in Fig. ^ We compare the frac¬ 
tional variability in our detrended lightcurves [aT) with the 
observed veiling fraction (Uobs). For A0620-00 and XN Vel 
1993 we use the detrended lightcurves with the relative phas¬ 
ing of the two sinusoids left free, but the difference between 
this and the fixed phasing model was negligible. For XN Mus 
1991, the ellipsoidal modulation is not adequately defined 
by the data so we do not detrend the lightcurves. Note that 
points from Zurita at al. (2002a), are typically for a redder 
bandpass than our V band measurements. For the veiling 
we use observed values based on spectroscopy, interpolating 
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Figure 7. a) The immediate field of XN Mus 1991, marked X, 
from one of our best images. The new variable is marked V. Stars 
Cl to C4 were combined as a reference for differential photometry 
of both stars. Stars C5 and C6 were used to check the extraction 
and estimate errors; C6 is shown on lightcurves as it is closer to 
the brightness of XN Mus 1991. b) Lightcurve of the new variable 
star. The same comparison star as was used for XN Mus 1991, C6, 
has been shown; this is somewhat brighter than the new variable. 


where more than one wavelength is available. The values we 
use are as collated by Zurita et al. (2002a), together with 
Orosz et al. (1996) and Casares et al. (1997) for XN Mus 
1991, and Filippenko et al. (1999) for XN Vel 1993. 

A precise correlation between veiling and fractional 
variability would obviously depend upon a number of as¬ 
sumptions and violation of these, together with the uncer¬ 
tainty in the measurements, will introduce scatter in the 
plot. The comparison can only be crude anyway, as the rms 
variability is not ideal for comparing datasets; it is effectively 
an integration of the PDS, but the limits of the integra- 
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Figure 8. Correlation between the fractional variability and the 
fractional contribution of disc light. Open circles are from Zu- 
rita et al. (2002a), filled circles are from this work. Triangles are 
spectroscopic observations from Hynes et al. (2002). The sources 
are: 1. A0620-00, 2. GS 2000+25, 3. V404 Cyg, 4. Cen X-4, 5. 
GRO J0422+32, 6. XN Mus 1991, 7. XN Vel 1993. The dashed 
line is a linear fit, passing through the origin, to all of the photo¬ 
metric points except that of XN Vel 1993. 
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tion depend on the length of the observation and the time- 
resolution, and hence will vary from dataset to dataset. The 
most serious physical assumption is that the veiling source 
is the same as the flaring source. This might not be the case, 
for example, if the veiling came from the whole disc, but the 
flaring only from the inner edge, or the stream-impact point. 
If the flaring source does differ from the veiling source, it 
could have a different spectrum, and hence it would be in¬ 
valid to combine R and V band observations in the same 
plot as we have done. There might also be differences in 
visibility, with, for example, the disc being foreshortened, 
but emission from an inner spherical flow being closer to 
Isotropic. This would lead to an inclination-dependent scat¬ 
ter. Finally intrinsic variations between sources will intro¬ 
duce further scatter. Consequently, it is not surprising that 
there is significant scatter. In spite of this, a correlation is 
visible; certainly for objects where the veiling (i.e. the disc 
contribution) is small, the variability amplitude is always 
small as expected. 

It is interesting to further extend the comparison be¬ 
yond optical photometry. Time resolved observations have 
also been performed in Ha for V404 Cyg (Hynes et al. 2002). 
These observations are useful because unlike the optical 
photometry, they should not show a significant contribu¬ 
tion from the companion star; the fractional disc contribu¬ 
tion is near unity. We have therefore added the fractional 
variability from these data to Fig. ^ for comparison. These 
points do lie relatively close to the line; they certainly con¬ 
tinue the sense of the extrapolation. This is not surprising 
given that Hynes et al. (2002) demonstrated that line and 
continuum flares in V404 Cyg are correlated. It Is possible 
that all optical variability in quiescent SXTs, continuum and 
line, has a relatively uniform undiluted fractional variabil¬ 
ity of 20-30 percent. In fact, the Chandra X-ray observation 


of V404 Cyg also showed a comparable level of variability 
(Kong et al. 2002). 

XN Vel 1993 does appear to show much less variability 
than expected from the very large veiling (60-70 percent) es¬ 
timated by Filippenko et al. (1999). Furthermore, their esti¬ 
mate was done at 6300 A, and we might expect a higher disc 
fraction in the V band. Such a large value does, however, 
seem difficult to reconcile with the pronounced ellipsoidal 
modulation which we observe. It could be that the disc was 
fainter at the time of our observations than when Filippenko 
et al. (1999) observed it, and hence that the veiling is less. 
We should also consider, however, that the veiling estimate 
for XN Vel 1993 is probably the least certain of those in 
the sample. It is faint, and the existing spectra were heavily 
blended with the bright star (more so than for our photom¬ 
etry). The spectral type is not well determined; Filippenko 
et al. (1999) use an MO template to estimate the veiling, 
but note that the spectral type could be as early as K6. The 
heavy blending may also result in some contamination of the 
spectrum of XN Vel 1993 by the brighter star, reducing the 
validity of the veiling determination further. Consequently, 
we cannot be confident that XN Vel 1993 does show signifi¬ 
cantly less variability than the trend suggested by the other 
sources in Fig. 


6.2 Lightcurves compared 

Since our Gemini observations of A0620-00 and XN Mus 
1991 show levels of variability consistent with their veiling, 
and they have similar orbital periods (7.8 and 10.4 hrs re¬ 
spectively), we might expect that the characteristics of the 
variability, such as the range of timescales present, should 
be similar; the only major differences between the observed 
variability properties should be in the amount of dilution 
by the non-varying light from the companion star. Some 
difference is already suggested by the PDS shown in Fig. ^ 
below 1 mHz, the PDS of A0620-00 appears to flatten, while 
that of XN Mus 1991 does not. Above ImHz, however, the 
PDS are similar apart from the differing normalisations. The 
PDS is never the whole story, however, as phase information 
is not preserved. We therefore should also directly compare 
the lightcurves. Superficially, at least, these look very dif¬ 
ferent. We must be a little careful, however; in XN Mus 
1991, the ellipsoidal modulation is weaker and the flaring 
much stronger than in A0620-00. The prolonged flares seen 
in XN Mus 1991 might, if scaled down in amplitude and 
superposed on a strong ellipsoidal modulation, appear to be 
undetectable or look like residuals to the fit to the modula¬ 
tion rather than real flares. Equally, the short, sharp flares 
in A0620-00 would be less striking when observed at lower 
time resolution. A more realistic comparison is therefore to 
rebin the A0620-00 lightcurve by a factor of 4 (giving a cy¬ 
cle time of about 32s, the same as used for XN Mus 1991), 
and scale the variations seen in XN Mus 1991 down to the 
same red-noise power and superpose them upon an ellip¬ 
soidal modulation. 

The results are shown in Fig. ^ After applying this 
processing to the lightcurve of XN Mus 1991, it does look 
more similar in character to that of A0620-00. The drop at 
the beginning of the lightcurve now blends indistinguishably 
with the synthetic ellipsoidal modulation. The strong, broad 
feature near the end is still visible, though less dramatic. At 
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Figure 9. A less biased comparison between the lightcurves of 
A0620-00 and XN Mus 1991. The lightcurve of A0620-00 has 
been rebinned to the same time resolution as that of XN Mus 
1991. The XN Mus 1991 lightcurve has been scaled down to what 
would be expected for the lower veiling of A0620—00, and an el¬ 
lipsoidal model appropriate to A0620—00 has been added. 

higher frequencies, however, one could easily believe that the 
two lightcurves are of the same source, in the same state. 
This analysis therefore supports the assertion that some of 
the very obvious differences between the lightcurves, at least 
of sources of similar orbital period, do arise from the differ¬ 
ing amplitudes of the ellipsoidal and flaring components, 
which in turn are simply dependent on the viewing geome¬ 
try. Differences in the low frequency PDS do suggest some 
real differences in the variability properties, however, and 
may give rise to some of the vertical scatter in the correla¬ 
tion plot (Fig. This does not necessarily require that the 
origins of the variability in A0620-00 and XN Mus 1991 dif¬ 
fer; both PDS can be accommodated in a broken power-law 
model, with the break out of the observable range in XN 
Mus 1991. 


6.3 Is the variability the same as that in 
cataclysmic variables? 

Further clues to the origin of the variability may be obtained 
through a comparison with cataclysmic variables (CVs). We 
might expect the outer regions of the disc to be rather sim¬ 
ilar in quiescent SXTs and in quiescent dwarf novae (DNe). 
Hence if the variability originates from magnetic reconnec¬ 
tion in the outer disc, or from the hot spot, the variabil¬ 
ity properties should be similar. Bruch (1992) has compiled 
rapid photometry of many CVs. If we select from his sam¬ 
ple only quiescent dwarf novae, a range of PDS slopes of 
— 1.6 to —2.6 is seen; the best studied case is SS Cyg which 
spans —2.0 to —2.6. These are systematically steeper than 
we see, but not by so much as to be conclusive; given only 
two well determined SXT PDS, we cannot rule out the possi¬ 
bility that these represent the flattest examples drawn from 
a similar distribution. The variations in quiescent DNe have 
B band full amplitudes of 0.18-1.26 mag. The latter is not 


directly comparable to the rms we measure, but for a Gaus¬ 
sian distribution of magnitudes would correspond to rms 
variations of 6-50 percent; this is a crude comparison, but 
at least enough to see that the amplitudes observed are 
similar to those in SXTs (Fig. In DNe, as in quiescent 
SXTs, a range of amplitudes is expected due to dilution 
of the variability by non-varying light, although unlike the 
SXT case, this usually comes from the white dwarf, not 
from the companion star. Unfortunately, existing observa¬ 
tions of quiescent DNe such as those of Bruch (1992) do 
not adequately sample low-frequencies to test for a break in 
the PDS, although a break has been seen in the VY Scl 
type CV, KR Aur during a high state (Kato, Ishioka & 
Uemura 2002). Thus the PDS of quiescent DNe and SXTs 
appear broadly similar, although DNe may exhibit a some¬ 
what steeper slope. This comparison is, however, inconclu¬ 
sive without a larger sample of well determined PDS from 
both classes of objects. 

6.4 The significance of band-limited noise 

The presence of band-limited noise in A0620-00, if not an 
artifact, is intriguing. This form of noise is also seen in 
low/hard state SXTs (see Wijnands & van der Klis 1999 
and references therein). This is illustrated in Fig. by 
comparison with the X-ray PDS of the low/hard state SXT 
XTE J1118-I-480 (Hynes et al. in preparation). The similar¬ 
ity of the PDS suggests that the optical variability could 
have a similar origin and might be associated with the cen¬ 
tral X-ray source, although we cannot say whether it rep¬ 
resents direct emission from an advective flow (e.g. self- 
absorbed synchrotron as predicted by Narayan et al. 1996) 
or results from heating of the outer disc. Like the quiescent 
state, the low/hard state is often interpreted as involving 
a truncated accretion disc with an evaporated central re¬ 
gion. An advective model has been successfully applied to 
XTE J1118-I-480 (Esin et al. 2001), and with spectral cov¬ 
erage into the EUV it was clear in this case that the disc 
must be truncated (Hynes et al. 2000; McClintock et al. 
2001). The most thorough analysis of this dataset indicates 
an inner disc radius of ~ 350Rsch (Chaty et al. 2002), sig¬ 
nificantly less than the 10^-10® Rsch usually invoked for the 
quiescent state (e.g. Narayan, Barret & McClintock 1997). 

There are differences in the two PDS shown in Fig. 

The normalisation is clearly lower in the case of A0620- 
00, as a fraction of the mean flux, but this is easily under¬ 
stood as its optical light includes a considerable contribution 
both from the companion star and other parts of the accre¬ 
tion flow which exhibit only very low frequency variability. 
There does not appear to be an analogue of the QPO seen 
in XTE J1118-I-480, but this is not always present in the 
low/hard state, either. The most significant difference ap¬ 
pears to be in the break frequency, which is much lower in 
A0620-00 than in XTE J1118-I-480. The origin of the break 
is not known, but it is plausible to expect that it approxi¬ 
mately scales with the size of the inner region. A low break 
frequency in quiescence would then be expected, since the 
advective region is expected to be larger. This trend is sup¬ 
ported by the much lower break frequency of 4 x 10~^ Hz 
seen in the AGN NGC 3516 (Edelson & Nandra 1999). The 
simplest expectation would be that the break frequency will 
scale linearly with the inner disc radius, but this may not 
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Figure 10. PDS of A0620-00 compared with XTE J1118+480 
in outburst (adapted from Hynes et al. in preparation). 


actually be the case; it will probably vary with some charac¬ 
teristic length scale, which could depend on the mass trans¬ 
fer rate or other parameters, as well as the inner disc ra¬ 
dius. Consequently there may be additional factors of order 
unity. If we assume that the relationship is linear, and that 
the black hole masses in A0620-00 and XTE J1118-1-480 are 
similar, then the ratio of break frequencies implies an inner 
disc radius in A0620-00 of ^ 1.1 x lO^^Rsch- For compari¬ 
son, advective models of quiescent SXTs, and specihcally of 
A0620-00, assume an inner radius of 10^-10“* Rsch (Narayan 
et al. 1997). Our estimate is obviously very crude, and could 
be off by a factor of a few, so the agreement with theoretical 
assumptions is reasonable. 

The absence of a break in XN Mus 1991, or its lower fre¬ 
quency, however, challenges this interpretation. If no break 
is seen to a factor of three lower in frequency then one 
would naively expect an inner disc radius three times larger, 
3 X 10“* Rsch- This seems large, at least for a relatively short 
period system such as XN Mus 1991. Consequently, if this 
interpretation of the break is correct, it seems likely that 
the break frequency does not scale linearly with radius, and 
that other factors do come into play. The absence of a clear 
break in XN Mus 1991 would make sense if the transition 
radius in A0620-00 was unusually small at the time of ob¬ 
servations. Indeed, the absence of a similar break in the PDS 
of Zurita et al. (2002a) requires some variation in the PDS. 
These considerations obviously reduce the potential value of 
the break frequency for measuring the transition radius, at 
least until the mechanism by which it arises is understood. 

One mechanism proposed for explaining the flaring in 
low/hard state systems involves a cellular automaton model 
for an accretion disc (or advective flow) in a self-organised 
critical state (Mineshige, Ouchi & Nishimori 1994a; Mi- 
neshige, Takeuchi & Nishimori 1994b). It is assumed that 
material is injected into a region subject to an instability if 
the density rises above a critical value; when this occurs an 
avalanche is triggered and the energy release is manifest as a 
flare. This model can reproduce a band-limited noise PDS. 


The break frequency is related to the size of the unstable 
part of the disc, which has subsequently been identified with 
an advective region (Takeuchi & Mineshige 1997). Assum¬ 
ing a 10 Modot black hole, a break frequency of ~ 10“® Hz 
implies the size of the region is ~ 3000 Rsch (Mineshige et 
al. 1994a), comparable to that obtained by scaling relative 
to XTE J1118-1-480. This radius also depends on the tem¬ 
perature and viscosity of the region, so is quite uncertain. 
Other aspects of the flare behaviour are consistent with such 
a model, such as the distribution of flare amplitudes and du¬ 
rations (Zurita e t al. 2002a) and the roughly symmetric flare 
profiles (Section 3.1) as modelled by Manmoto et al. (1996). 


7 CONCLUSIONS 

We have studied short-timescale flaring in three quiescent 
SXTs. This flaring is detected in all three sources, at low 
amplitudes in A0620-00 and XN Vel 1993, and at a much 
higher level in XN Mus 1991. With the large aperture of 
Gemini it is possible to observe at higher time resolution 
than in previous studies. We find that the variability ex¬ 
tends to the shortest timescales observable, with pronounced 
changes sometimes seen in 30 s or less. This is also indicated 
by the extension of the red noise component in the PDS to 
0.05 Hz or even higher. 

A comparison of our observations with those of Zurita 
et al. (2002a) supports their conclusion that the flares are 
associated with the accretion flow rather than with the com¬ 
panion star. This is clearly shown by a correlation between 
the variability amplitude and the fractional disc contribu¬ 
tion to the spectrum. The amplitude of Ha variations in 
V404 Cyg is also consistent with the correlation, and it is 
likely that the line and continuum variations have a related 
origin. 

Compared to quiescent DNe, i.e. the nearest compara¬ 
ble systems not containing a black hole or neutron star, our 
sources show similar levels of variability, but perhaps with 
a somewhat flatter PDS. A rigorous comparison will, how¬ 
ever, require a larger sample of objects and more intensive 
observations. 

In A0620-00, we detect a low-frequency break in the 
PDS at ~ 10“® Hz. The PDS overall looks very similar to 
those of low/hard state SXTs. If the break frequency scales 
linearly with the size of an inner evaporated region, then 
a comparison with XTE J1118-1-480 suggests that this re¬ 
gion has size ~ 10^ Rsch in A0620-00 in quiescence, al¬ 
though there are likely to be other factors involved and 
this is an extremely crude estimate. No such break is con¬ 
fidently detected in XN Mus 1991; if present it must be at 
< 3 X10“^ Hz. It may be that the break frequency varies from 
source to source and epoch to epoch, possibly in response 
to changes in the inner truncation radius of the outer disc. 
Again, more intensive observations are needed to explore 
this behaviour. 
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